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Abstract 

To study the rate of energy and momentum loss of a heavy quark in QGP, specif- 
ically in the hydrodynamic regime, we use fluid/gravity duality and construct a 
perturbative procedure to find the string solution in gravity side. We show that by 
this construction, the position of the world-sheet horizon and thereby the drag force 
exerted on the quark can be computed perturbatively, order by order in a boundary 
derivative expansion, for a wide range of quark velocities. At ideal order, our result is 
just the localized expression of the drag force exerted on a moving quark in thermal 
plasma, while for a quark whose velocity does not belong to the mentioned range, we 
predict a nonlocal drag force. Furthermore, we apply this procedure to a transverse 
quark in Bjorken flow and compute the first-derivative corrections, namely viscous 
corrections, to the drag force. We show that the diffusion time depends on the initial 
momentum of the transverse quark and can be about 4fm for a high energy quark 
which has been created sufficiently late. 



*navid_abbasi@physics.sharif.ir, abbasi@ipni.ir 
^davody@ipm.ir 



Contents 



1 Introduction [I| 

2 World-sheet covariant description of drag force ^ 

3 Moving quark in a boosted thermal plasma |^ 

3.1 Drag force computation ^ 

3.2 Embedding of the string H 

3.3 Covariant results ^ 

4 Gravity dual to a quark in fluid P 

4.1 Review of fluid/gravity duality |1^ 



4.2 Construction of perturbative procedure |ll 



5 Quark in temperature varying plasma ^ 



6 Quark in Bjorken flow 18 



6.1 Perturbative computations 11^ 



6.2 Phenomenological estimates ^ 



7 Discussion B7 



A Appendices ^8 

A.l Global boosted plasma p 



A. 2 Fluid dynamical correction ED 



A.2.1 Free fluid 29 



A.2.2 forced fluid EDI 



A. 3 Varying plasma ^ 



A.4 Bjorken flow 30 



1 Introduction 

Studying systems out of equilibrium is an interesting subject in physics. Depending on 
how it evolves towards equilibrium, there are diflerent methods to study the system. In a 
familiar case, when the deviation from equilibrium has small amplitude everywhere, one 
can apply linear response theory to study the dynamics of the system. Computing retarded 
Green's functions in equilibrium, we can extract certain information about the dynamics of 
the system slightly deviated from its equilibrium state. However, when the understudying 
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system is described by a strongly coupled field theory, we are not able to do such compu- 
tations through perturbative methods. In such cases, AdS/CFT correspondence appears 
powerful tool. 

AdS/CFT conjecture in its original statement, relates the type IIB string theory on 
AdSs X space-time to the four- dimensional M = 4 SYM gauge theory [Q. On the other 
hand, there are many examples which have been studied with the holographic description 
of AdS/CFT in which, a strongly coupled field theory living on the boundary of the AdS 
space is pictured to the weakly coupled gravity theory in the bulk of AdS . Accordingly, 
this holographic picture, provides a tool to study the strongly coupled systems out of 
equilibrium 

Although our familiar strongly coupled systems do not exactly coincide with the A/" = 4 
SYM gauge theory, we apply this theory as a useful model to study such systems through 
AdS/CFT correspondence. This modelling shows agreement with real physical systems in 
particular circumstances. For example in QGP experiments at RHIC or at LHC, a gas 
of quarks and gluons is produced with the temperature about 170 MeV which is strongly 
coupled. Although QCD as the underlying theory in these experiment is distinct from 
A/" = 4 SYM gauge theory, but this has not prevented us to apply this to to study of 
QCD in QGP experiments. Specially in the hydrodynamic regime, computation of the 
transport coefficients is known as a remarkable success of AdS/CFT providing results 
in agreement with experiment. 

One of the other significant subjects from both experimental and theoretical stand- 
points, is the motion of quarks created during heavy-ion collisions in the plasma. Up to 
now, invoking AdS/CFT, three different mechanisms have been introduced to describe the 
rate of energy and momentum loss through this motion; exerting drag force on the quark 
from the plasma |^, Cherenkov radiation of a fast heavy quark and radiation related 
to quark's decleration P|. There is a simple holographic picture for these problems. The 
probe quark moving in plasma is mapped to a probe string in the AdS space. So instead 
of taking into account the quark motion in a strongly coupled system one can simply trace 
the classical dynamics of a string in the gravity side. We will concentrate on the first case 
in this paper. A common part in both related seminal works P, ^ is that the authors have 
studied a the uniform motion of a quark in the boundary plasma. This quark is the end 
point of a string stretching in the AdS bulk. So in order to move uniformly, it should be 
forced by an external source, i.e. an electric field. The momentum rate fiowing down to 
the string has been interpreted as drag force exerted from the plasma on the quark. It has 
been shown in ^ that when a quark is dragged with a constant velocity v in a. thermal 
plasma at temperature T, the rate of the energy and momentum loss of the quark in the 
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plasma frame are given by 



dp 1 dE IT ^r, V , , 

= - — = 1.1 

dt V dt 2a' VI 

In these derivations, besed on AdS/CFT correspondence, the gravity set up dual to the 
thermal field theory has been taken as a uniform 5-dimensional AdS black brane with the 
metric 

dT*^ 

ds^ = GM^dx^dx^ = — — - + r\P^, - f{hr)u^,u,)dxW (1.2) 

where M = {r,fi), f{r) = 1 — ^ and the temperature of the black brane is given by 
T = l/nh. Also P^u = u^Uy + Tj^y is a projection tensor where is the four- velocity of 
the plasma. For a plasma at rest in the laboratory frame = (1, 0, 0, 0). 

From the experimental viewpoint, an important case to study is a moving quark in an 
expanding plasma where, local thermal equilibrium has been attained. It can be generated 
in the time interval between thermalization and hadronization, after scattering of two 



massive nuclei [|^. In this time interval, the excited matter is expanding and cooling 
down such that at any point, hydrodynamic description is valid and therefore, dynamics of 
the expansion can be described by the corresponding equations of motion [|lT]. Beforehand 
AdS/CFT, there was not any theoretical approach to study the hydrodynamic regime in 
strongly coupled field theories. In fact, AdS/CFT has opened a new window to research: 



study of CFT hydrodynamic from gravity. First efforts in this context were made by |T^, |T3 
which were confirmation of existence of the hydrodynamic regime in a A/" = 4 SYM gauge 
theory. In another remarkable development, Janik and Peschanski found an asymptotic 
gravity dual to the hydrodynamical fiow in QGP, namely Bjorken fiow, from a nonlinear 
derivation \T^\. In 2007, BHMR introduced the full version of fiuid/gravity duality [T3 



in which, they showed that there is a one to one map between fiuid dynamical fiows of a 
CFT on the boundary with the long wavelength perturbations of the AdS black brane in 
the bulk. 

It is natural to investigate the quark motion in a dynamical fiow as a phenomenological 
problem related to the QGP experiments. In this direction, the drag force exerted on a 



quark moving in Bjorken fiow has been computed in and later in ^^7j. In these works, 
Bjorken fiuid has been considered as an ideal fiuid. To be more precise, one can go further 
in the hydrodynamic expansion and study the viscous effects. For example, by applying 
the gravity set up found in [0 , it is possible to compute the drag force exerted on a quark 
moving in viscous Bjorken fluid. 

As discussed, fluid/gravity correspondence maps each long wavelength perturbation 
around an AdS black brane to a hydrodynamical flow on the boundary of AdS. Apart 
from sketching the mentioned map, a key point in fluid/gravity duality is that it provides 
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a perturbative technology to find the regular solutions in AdS space. Similarly, based 
on fluid/gravity duality, we construct a perturbative procedure to find the classical string 
solution in the long wavelength perturbed AdS background. The choice of coordinates has a 
key role in constructing our perturbative procedure. We exploit the Eddington-Finkelstein 
coordinates in which the extension of the string in the boundary directions is finite. It is in 
contrast to the coordinates systems taken in |^ where, the embedding of the string has 
an infinite extension in the boundary spatial directions. This delicate point is the origin 
of our perturbative procedure. 

It is important to note that constructing of such perurbative procedure is the key issue 
in development of our idea in this work. For this procedure to be applicable, it is necessary 
that the relative velocity of fiuid and quark to be lower than a relativistic upper-bound 
velocity which turns out to be about 0.97. Using this procedure, one can compute certain 
physical quantities related to a quark in fiuid medium perturbatively. For example we 
proceed to compute the drag force exerted on a quark moving in a fiuid dynamical fiow. 
As a main result, we find that at ideal order, it can be computed in a fiow- independent 
manner. In addition, we extend this procedure to higher orders of perturbation. We show 
that in a wide range of the quark velocities, one can compute corrections to the drag 
force perturbatively, order by order in a boundary derivative expansion. Finally, we apply 
this method to a transverse quark created in Bjorken fiuid and compute the drag force 
perturbatively up to viscous order. Our numerical computing shows that the diffusion 
time depends on the creation time of the quark, with respect to the collision time. For a 
transverse quark created at r = 15/m, we estimate it about 4/m. This is an improvement 
of the result of |jl9| where compared to the result of |^0|]. We have considered a varying 
temperature profile through the motion which was lacking in previous approaches. 

This paper is organized as follows. We begin the section 2 with a brief review of drag 
force derivation in String/ CFT duality in a world-sheet covariant approach. In § |] we apply 
the statement of § ^ to compute the space-time covariant drag force exerted on a quark in 
a global boosted plasma in. § § is devoted to the core of this paper in which after reviewing 
some key points in fiuid/gravity correspondence, we explain our perturbative procedure. 
Then, we apply this method to find classical string solution in gravity background dual 
to a boundary fiuid dynamical fiow. As a preliminary application of our perturbation 
method, we study a quark moving in a slowly time varying plasma in § ^. In § ^ we apply 
our perturbation to Bjorken fiuid and compute first-derivative correction to the drag force 
exerted on a transverse quark in this fiuid. In particular, we find the diffusion time of the 
quark. Finally in § ^ we end with a discussion of results and of possible future projects. 
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2 World-sheet covariant description of drag force 



Consider a quark moving with velocitywq in the boundary field theory. According to the 
string/CFT duahty, such quark is the end point of an open string, embedded in a one 
higher dimensional space-time dual to the boundary field theory. The dynamics of the 
classical string is described by the Nambu-Goto action as follows 

S = -^ jdadre^V^, ae[0,a,] (2.1) 

where $ is the background dilaton field. Also in the above expression, g = detgajs where 
gai3 = Gfiu daX^dpx" . Free motion of the quark in the plasma is dissipative. AdS/CFT 
duality ascribes this dissipation effect to flowing the momentum down to the string in the 
bulk. Therefore, to have a uniform motion, the quark should be forced by an external 
source such as an electric field on the boundary. In such a case, the mentioned fiow rate is 
interpreted as the drag force exerted on the quark. 

By construction, drag force should is a vector in space-time and a scalar on the world- 
sheet. In order to have it in a world-sheet covariant form, consider the world-sheet currents 
as following 

dC e5*(^'') 

where a and M = r, fi describe the world-sheet and space-time directions respectively. In 
contrast to P^, these currents are not the world-sheet vectors. According to the previous 
paragraph discussion, to compute drag force we should find ingoing fiux of the energy and 
momentum across the time-like boundary of the world-sheet, the a = ai curve. Assuming 
n° as the inward normal vector to this curve, the covariant 4-drag force is given by: 

F, = -P; n«. (2.3) 

It is clear that the choice of the gauge determines what combination of P^s is equal to the 
drag force. Let us recall that the relativistic four-force, = -^p^, exerted on the quark 
should be perpendicular to the quark four-velocity: 

F, < = 0. (2.4) 

For a quark moving with velocity Ug in the 1-direction of the plasma frame, applying p.4| ) 
to = 7 (^, f , 0, 0) ensures that 

? = (2.5) 
at Uq at 
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which can be obviously observed in ( |1 . 1| ) and all of our perturbative results in next sections. 



3 Moving quark in a boosted thermal plasma 



Transforming the coordinates in AdS black brane metric ( |1.2|) to Eddington-Finkelstein 
coordinates, we find the metric of gravity space-time dual to the boundary boosted thermal 
plasma: 

ds^ = GMNdx^'^dx^ = -2M^dx^dr + (P^^ - f{br)u^u^) dx^dx"" (3.1) 

where, the boosted plasma with velocity and temperature T = ^ lives on the boundary 
of space-time, r — )■ oo. It is significant to note that (|3.1D introduces an exact 4 + 1 



parameters family of Einstein equations solutions ^1 \ . 

To study the energy and momentum loss of a moving quark in this system, we consider 
the four-velocity of plasma to be = (m°, u^, 0, 0) and for simplicity assume that the quark 
is moving in the 1-direction with a constant velocity Uq. Above-mentioned quark is the 
end point of an open string embedded in a 5-dimensional space-time by the metric (|3.1|) . 
Dynamics of the string is governed by the Nambu-Goto action ( p.l|) where dilaton factor 
can be dropped in the above background. In the static gauge (a, r) = (r, t), embedding of 
the string is given by X^'\r,t) = {r,t, X^{r,t)). Plugging 

= Uqt + ^{r), X^ = X^ = (3.2) 

in (|2.1|), we reach to following expression: 



C = V'Af2T25fTC' (3.3) 

where A, B and C are some constant coefficients which depend on the bulk metric elements 
and the velocity of the quark ( |A.1|) . As it can be observed, ^ enters into the world-sheet 



lagrangian only through its derivative ^' . So instead of solving EOM, one can exploit the 
conservation equation of the conjugate momentum of ^: vr^ = dC/d^' . Deriving ^' in terms 
of vTg, we find 

B 1 / ^B^-AC , , 



It can be simply checked that B^ — AC vanishes at some point in the bulk. Obviously, for 
the second term to be real and so the induced metric to be non-degenerate, numerator and 
denominator should have common root, namely r*. Solving B^ — AC = 0, r* is obtained 
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as 

Demanding the denominator to vanish at r*, determines the value of vr^ as: 

where /3 = , is the 1-component of the plasma velocity in the boundary lab frame 



(LF). Inserting ( |3.6| ) in (p.4p, we have: 

^1Z^^(_1± 1 (3.7) 

r2/(6r)v/r^^ b^r^ 

where the choice of + or — depends on the sign of Ug — (3; for Ug > (3 {ug < f3), the + (— ) 
is acceptable. 

3.1 Drag force computation 

Let us recall from of our discussion in § ^ that the above gauge choice, fixes the inward 
normal vector to the boundary of the world-sheet at r — )■ oo as tt-q = (—1,0). So the 
1-component of 4-drag force may be written as: 



= Fi = PI = (^-^u{^ 



(3.1 



where 

1 ^ (x.x)(x7- (x)2(x'^)' 1 



= Gi / • \ ^ ' ^ = TT.. (3.9) 

^ 27ra' v/=^ 27ra' ^ ^ ^ 



Translational invariance of (|3.3|) ensures the conservation of the world-sheet currents, 
dall'^ = 0. Knowing the value of 11^ at r*, one can integrate this equation to compute 
as a function of radial coordinate. Since the space-time is static, it turns out that 11^ 
is constant along the string and equal to 11^ (r*). Therefore we obtain the 1-component of 
4-drag, F>^ = (7^,7^,0,0), as follows: 

. 1 dp 1 1 , ^ 

F = , -T- = n 3.10 

^/T^g dt 2'Ka' v^r^wf 

where dp/dt is the 1-component of 3-drag force in LF. 

It can be simply seen that for Ug > P, drag force is pointed to opposite direction of the 
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quark motion while for Uq < (3 aligns in the quark direction of motion. Expectedly, when 
the quark is moving with the velocity equal to that of the plasma, does not sense any drag 
force. 

To compute the rate of energy loss, we should work out Tl\. After a simple calculation 
we have 

= A— ¥ = r^^,. (3.11) 

It is interesting to note that for /3 = 0, the rate of energy and momentum loss in the LF 
obtained from ( |3.10D and ( ^.11| ) exactly coincide with ( |1 . 1\) . It is expectable that one could 



obtain the above expression by a lorentz transformation of (|1 . 1|) . In addition, ( p. 101 ) and 
( p.liP remind us that instead of studying string picture dual to a moving quark in global 
thermal plasma, one can apply string/CFT duality to a rest quark in a global boosted 
thermal plasma, since both {uq = v, P = 0) and {uq = 0, (3 = —v) cases have common 
results from (13. ( 



3.2 Embedding of the string 

The embedding of the string in the bulk space-time is the solution of equation 
Assuming Uq > P and picking up the plus sign to ensure the momentum flowing down the 
string, we have 

e(r) = -^7== b {tan-\rb) - ^) (3.12) 

where we have chosen the constant of integration such that to be identified hj x = vt 
on the boundary. In Fig.|l| the embedding of the string has been compared in two different 
coordinates systems. Figure (a) shows the profile of the string for Uq > (3 in a. constant time 
slice of the metric ( p..2| ). In case (b) we have shown the same string has been considered 
in the coordinates system given in (|3.1| ). Although the quark is moving to the right, the 
string trails out in front of the quark in this case. It is not so surprising because, the time 
coordinate in ( p.l| ) is the ingoing Eddington-Finkelstein time in which all of events on an 
ingoing null geodesic are simultaneous. 

In contrast to the former case in which the extension of the string in the 1-direction 
is infinite, in the latter case, the string has occupied a finite range in the 1-direction. 
Substituting the value of the horizon radius in ( p.l2| ) gives the length of this finite extension 
in the LF as following: 

Ax^ = -b — ^==. 3.13 

The finiteness of Ax^, as a direct consequence of applying the Eddington-Finkelstein 
coordinates will play a key role in study of moving quark in fluid medium in next sections. 
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Figure 1: Comparison of the string embedding in the bulk between two different coordinates 
systems. In the Eddington-Finkelstein coordinates, string has been stretched in front of 
the moving quark in a finite range on the boundary. 

3.3 Covariant results 

In this subsection we want to present our previous result in the frame-independent expres- 
sions. It should be noticed that we have only two independent four- vectors in our analysis, 
four-velocity of the quark and of the plasma. So to have a covariant drag force we con- 
struct a linear combination of these velocities and demand the perpendicularity condition 
in (U). We obtain 



= 7g(^' 0' 0)' applying the LF results given in (pTTI). 
The position of the world-sheet horizon, r*, as a lorentz scalar can be given in following 
covariant expression: 



4 Gravity dual to a quark in fluid 

Our final goal in two next sections is to investigate the uniform motion of a quark in 
a conformal fluid described by the hydrodynamic regime of an interacting field theory. 
According to the string/CFT picture, to study the quark motion in such an interacting 
field theory, the string should be solved in a bulk gravity background dual to the boundary 
fluid flow with appropriate boundary conditions. So it is needed to recall some of the main 
results of fluid/gravity correspondence. 




(3.14) 





(3.15) 
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4.1 Review of fluid/gravity duality 

Consider a 4- dimensional field theory living on the boundary of an asymptotically AdS^ 
space-time. It has been shown that there is a map between the hydrodynamic regime 
of the strongly coupled boundary field theory and a class of inhomogeneous, dynamical 



black hole solutions in AdS space-time ^ . In fact the long wavelength perturbations 



of an asymptotically AdS black brane are mapped to corresponding dynamical flows of a 
conformal fluid on the boundary. 



According to |15], the dual gravity of an arbitrary boundary flow can be found pertur- 
bativly, order by order in a boundary derivative expansion. It means that the gravity dual 
to a boundary fluid flow with velocity u^{x°') and temperature T(x'^) can be obtained per- 
turbatively, order by order by adding appropriate correction terms to the localized version 
of the metric (B.ll) : 



ds^ = -2n^(x")dx^dr + r^[P^^{x'^) - /(6(x")r)u^(x")u^(x")]dx'^dx'^ + ds^, (4.1) 
where P^i^ = r]^y + u^u^ is the projection tensor. The first-order correction of the metric 



has been given in |A.2.1 



Although the above gravity set up seems so general, one may be interested to take into 
account forced fluid flows rather than such free flows on the boundary. This is the subject 
of in which, the bulk solution dual to an arbitrary fluid dynamical flow of the boundary 
field theory with an arbitrary slowly varying coupling, living on a weakly curved manifold, 
has been found up to second order in boundary derivative expansion. In order to describe 
such dynamical flows, instead of S* = j C the action of the boundary field theory is in the 
form S = J y/ge~'^C where gfj_i, is the weakly curved metric of the boundary and (f){x^) is 
an arbitrary slowly varying function. From this action, the Navier-Stokes equations will 
appear with a forcing term in the right hand side: 

V^T'^^ = e""^ C V>. (4.2) 



In , it has been explained that the gravity dual to an arbitrary general fluid dynamical 
flow in this field theory is a long wavelength solution of the Einstein-dilaton system with 
appropriate boundary conditions. Note that for dilaton field $(r, x'^) in the bulk, the 
boundary condition is the necessity of this field to be asymptoted to the given slowly 
varying function (j){x'^) on the boundary. Requiring expected boundary conditions, both 
metric of the bulk and ^{x'^) can be computed perturbatively, order by order in boundary 
derivative expansions. This perturbative procedure adds derivative correction terms to 
zero-order expressions of the metric and the dialton field in each order of perturbation: 

ds^ = -2u^(x")dx'^dr + [P^^(x") - /(6(x")r)u^(x")u^(x")]dx^dx'^ + dsl 
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$(r, x") = (pix") + $eor(r, x") (4.3) 

where Pfj_„ = g^y + u^Uy is the projection tensor. The first-order correction of the metric 
and the dilaton have been given in |A.2.2| . 



Emphasised by the authors of |15|, a key point in all above derivations is the appli- 
cation of ingoing Eddington-Finkelstein coordinates which in addition to a transparent 
presentation of horizon regularity, provide a clean picture of bulk gravity dual to the lo- 
cally equilibrated fluid domains on the boundary. In this picture, each locally equilibrated 
domain in the boundary fluid, extends along the ingoing null geodesies in an entire tube 
into the bulk. The bulk space-time within each tube is approximately a uniform AdS black 
brane where the width size of these tubes is the scale of variations in the fluid, namely L. 
As we know, having a gravity perturbative solution in boundary derivative expansion is 
possible, provided 

LT > 1 (4.4) 

at each point on the boundary. T is the local temperature at the point. 

Fluid/gravity duality has a speciflc characteristic. In each order of perturbation the 
relevant gravity solution is dual to the fluid dynamical flow in one lower order. Nevertheless, 
the stress tensor of the fluid will be determined at the same order to which the metric has 
been solved. Let us give an example. The flrst-order gravity solution is dual to an ideal 
fluid flow (zero order in hydrodynamic expansion) on the boundary but it provides the 
viscous fluid energy momentum tensor (flrst-order in hydrodynamic expansion). ||15| . 



4.2 Construction of perturbative procedure 

After the above brief review of fluid/gravity duality, we are going to study the quark mo- 
tion in a fluid background. Going from global equilibration state of the plasma to the 
hydrodynamic regime, the physical quantities related to the quark motion get hydrody- 
namic corrections. Meanwhile, the dual string picture of the quark changes accordingly. 
The main idea of our paper is that for a wide range of quark velocities, one can compute 
these corrections perturbatively, order by order in a boundary derivative expansion. 

The key point that allow us to use hydrodynamic perturbation procedure is exploiting 
the Eddington-Finkelstein coordinates by which, the whole string extension in the bulk 
may lie within just one tube. If such an extension exists, we can use the ultralocal method 
described in [|l^]to flnd the string solution perturbatively. To do that, we start with a string 
extended in the bulk and restricted to a tube with size 1/T. At zero order of perturbation, 
this assumption allows us to neglect the velocity and temperature variations on the the 
boundary slice of the tube, namely the boundary patch. Simply speaking, the difference in 
values of hydrodynamic quantities between two different points in a patch, is in the next 
order of perturbation. 
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So far, we have described a string embedded in a uniform AdS black brane space-time 
and therefore the drag force can be computed just hke as was done in § ^ where, the velocity 
of plasma can be chosen equal to the fluid velocity in an arbitrary point in the boundary 
patch i.e. in the position of the quark. So, as far as we can neglect the variation of the 
hydrodynamic quantities in a patch, the drag force exerted on a quark in a fluid can be 
computed by limiting (|3.14|) where the fluid global quantities are replaced by the local ones 



in the quark position. Consequently, we have 

F'i^n = i ^ [{^.-M^^K + ^'fi^n] (4.5) 

where is the fluid four-velocity and x" are the coordinates of the quark position. Since 
the local profile of the fluid has entered, the above expression identifies the drag force at 
ideal order. It is exactly similar to identification of the ideal fluid stress tensor in |[15|. 



Computing the drag force at the level of viscous fluid needs to considering derivative 
corrections. Exploiting our perturbative procedure, we proceed to find the string solution 
in first-order corrected gravity background as the dual of a quark motion in a viscous fluid. 
It should be noted that our constraint of the string stay in a single tube is highly restrictive. 
This strict condition emphasises that to compute the hydrodynamic corrections to string 
embedding, one should choose the rest frame of the quark (RF) where this restriction has 
a manifest physical meaning and can be easily implemented in. 

Now, we want to explore that under which circumstances in boundary side, the whole 
string extension stays in a single tube? Consider a one directional fluid dynamical flow 
in the LF. At zero order, where the velocity and temperature are global in a patch, the 
extension of the string in 1-direction of the boundary is given by ( |3.13| ) and in the RF in 
the fluid we have 

Ax' = ^ bix-) u\x-) = - ^f'K) (4.6) 

where u' = ^j^-^ position of the quark. 

It is clear from the above discussion that as far as the string extension is of order 1/T, 
the string lies in a tube. However ( [4. 61) shows that the restriction of the string to stay in 
a single tube, limits the relative quark-fluid velocity to values lower than an upper bound. 
Demanding | Ax^ | < 1/T, one can find 

/3 < ~ 0.97. (4.7) 

x/T7 



This upper bound for the fluid velocity in translates to the upper bound of quark velocity 
in LAB 

Vu.h. ~ 0.97 (4.8) 
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Figure 2: String in a bulk tube. In the RF and in the Eddington-finkelstein coordinates, 
the trail of the string is in opposite direction of the flow. 



and Since it turns out to be very close to the velocity of light, one need not worry about 
applying our perturbative method to relativistic cases. In practical cases, the charm quarks 
created with energies less than 6Gev satisfy our criteria. 

As a result, we have shown that when the relative velocity of quark and fluid is lower 
than ( [4. 81) , drag force can be expressed by ( [4.5|) at ideal order. This expression is indepen- 
dent of fluid dynamical flow and can be applied to any arbitrary flow which is compatible 
with fluid dynamics equations of motion. This result is in perfect agreement with |16, p!7| , 
where the computed drag force exerted on an external quark moving in the ideal Bjorken 
fluid is equal to the localized version of the drag force in global equilibrated plasma. But 
our construction shows that there is a concrete restriction on applying such local drag 
force: being the velocity of quark smaller than the upper bound (|4.8|) . 

What can we say about a fast quark with a velocity bigger than the above relativistic 
bound? In this case, the different points on the string do not lie in a boundary patch. 
Therefore, the drag force can not be a local function of the quark position which means that 
the drag force at some point in the fluid is a function of fluid dynamical fields (temperature 
and velocity) at some other point. This nonlocality is an AdS/CFT prediction which has 
not been obtained yet in the field theory side. As in pointed, the part of string above 
the world sheet horizon in the bulk corresponds to the "gluon cloud" of quark in plasma. 
Going with a velocity bigger than the upper bound leads to exiting the gluon cloud from 
the boundary fluid patch. Quark's gluon cloud lags behind it and so the nonlocal drag 
force appears. This is just an intuitive interpretation of this effect in the filed theory side, 
this point requires further investigation which we leave to future works. 

To compute viscous corrections to drag force, we make attempt to find the string 
solution in the first-order corrected gravity by applying our perturbative procedure. It will 
be the subject of the two next sections. 
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5 Quark in temperature varying plasma 



In order to compute derivative corrections to drag force exerted on a quark in fluid medium, 
we first choose a simple but instructive example. We will see that this toy model problem 
gives more insight into our perturbative procedure and reveals some of the ingredients 



discussed in 5 4.2 



Let us consider a quark uniformly moving in the 1-direction of a globally equilibrated 
plasma. If the temperature in the LF is globally varying in time, it seems to be as a 
dynamical flow with the following profile: 

= (1,0,0,0), T{t). (5.1) 

But as it can be simply seen, this profile does not solve the free fluid equations of motion, 
V fjT^'^ = 0. Instead, applying an appropriate varying dilaton function, this profile can 
be accounted as a forced fluid flow on a flat background. According to [^, the forced 



Navier-Stokes equations, ( |4.2| ), at the level of leading contribution of the dilaton function 
is given by: 

V^T^^ = -i7TT)^d''(l)do<P (5.2) 

where, at ideal order, the dilaton field does not give any contribution and therefore the 
temperature remains constant in time. So, the given profile in ( p.l|) changes to: 



= (1,0,0,0), T = Const, 0(t) (5.3) 

which says that at ideal order, a rest fluid with varying temperature does not exist. Instead, 
we could have a thermal field theory with a slowly time varying dilaton function. From the 
previous section, we recall that the drag force at ideal order can be obtained by localizing 
the thermodynamic drag force. Restoring dilaton factor in the computations of § p.l| , we 
can rewrite the rate of momentum loss in the fluid frame as 

1 ^ dp ^ ^ .^(,) ^2 (5 4) 

^ dt 2a' v^T^^ 



where from [23|, ^(t) is equal to (j){t) at zero order of gravity set up. From the string 
theory we remember that in the presence of a nontrivial dilaton field, the string coupling 
would be changed by a dilaton factor and so the 'tHooft coupling |^ 

A = glj^ N, = ATig, N, = \ (5.5) 
would be changed too. We see that requesting to have a varying temperature plasma is 
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equivalent to varying the coupling in this medium. Consequently equation ( |5.4| ) gives drag 
force exerted on a moving quark in a plasma whose coupling is slowly varying in time. 

Now, let us consider a quark in this medium and investigate the first-derivative cor- 
rections, namely viscous corrections, to drag force. According to the previous section 
discussion, an appropriate frame to study the derivative corrections is the RF, so hereafter 
we assume the slowly time variation of dilaton to be occurred at this frame. At viscous 
order, the profile of such flow is given by 

M^ = (n°,n\0,0), T(r), </>(r) (5.6) 



where = — 7m„ and r is the time coordinate in the RF. Based on [|23[, the dual picture of 



such flow in gravity side is an Einstein-dilaton system whose metric does not get correction 
up to the first order and first-order correction to the dilaton field is given by 



Jrb 



As before in the static gauge (o", r) = (r, r), the embedding of the string may be written 
as X^{r, t) = (r, r, X\r, r)). Taking the following ansatz 

X^ = ^(r,r), X^ = X^ = 0, (5.8) 



the Nambu-Goto lagrangian will be given by 

C = e^*("'") ^A^'^ + 2B^' + C. (5.9) 

In the above expression, A = Ai, B = Bi + B2C,, C = Ci + 2C2C, + Cs^"^, where among the 
subscripted coefficients only Ai is a function of time and the others are constant ( |A.3| ). 
Also ^' and ^ are partial derivatives of ^ with respect to r and r. 

As discussed in § ^TT] we should find the world-sheet currents to compute the drag force. 
Based on our discussion in § ^ drag force in the present case is given by ( ^.8|) similarly. 
From (|5^ ), it's clear that H^s are not all conserved and we have 

a„n^_2,3 = 0, d^w^ = u. (5.10) 

Note that a = r, r are the world-sheet coordinates and yU = r, 1,2,3 are the boundary 
directions. The existence of time dependency in the action, does not allow 11^ to be 
constant along the string. Physically, this variation along the string is related to the 
energy and momentum transfer between the string and the non-static bulk of the space- 
time. Thus by integrating the equation ^qII" = over the r component of the world-sheet 
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we obtain: 



Ul{r,T)-U[{ro,T) = -A.I\l{r',T)dr' [5.11] 



where 



ni(r, r) = --i^ei^(-^-)G,. (''-^^(^'^^;^(^^'(''^^^ = --i^vr(r, r). (5.12) 
Zna \ —Q ZTxa 



To study the string embedding, we derive ^' from the above equahty 



(5,13) 

with y4, i? and C introduced in ( |5.9|) . Ahhough the presence of a slowly varying dilaton 
function have changed the numerator and denominator in comparison with ( |3.4|) , just like 
there we demand the second term in (|5.13|) to be real . Again, this reality condition obliges 
the numerator and denominator to have common root. To compute this root perturbatively, 
we add a correction term to the embedding as 

e(r,r)=eo(r)+U(r,r). (5.14) 

^o('") is the zero-order embedding given in ( ^.12p in which Uq vanishes in RF: 



Wr)=u)h (tan"i(r&) - ^) . (5.15) 

As a remarkable point, the above expression is independent of $; something that can 
be simply traced in § ^. The other point is the constancy of h at this order. It is a 
direct consequence of equation (p.2|). Inserting ( |5.14| ) in ( |5.13| ) and Solving — AC = 







up to first order, we obtain the first-order corrected position of the world-sheet horizon 
\\f^f ^^^^ presence of such a dilaton function has not changed 

r* up to the first order. Confirming our expectation that a scalar field does not influences 
the geometry at the order. In order for r* to be also the root of the denominator up to the 
first order, we should have: 



oo ^3 1 \ „i (5-16) 



As clearly observed, the only viscous contribution to above expression is related to the 
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first-order correction of $ given in ( [5.7| ). 

Up to now, we liave succeeded to compute tlie value of 11^ (r,r) at a special point on 
the string. In fact this point is the first-order corrected horizon of the world-sheet which 
helps us to compute 11^(0, r) |r^oo from ( [5.11D . To compute the right hand side of ( [5.11D 
up to the first order, the integrand should be worked out at zero order: 



27ra' 



1 + rW 



Inserting ( |5.16|) and (|5.17|) in with choosing tq = r*, we obtain 11^ (r,r) up to the 

first order as following: 



n^(r,r) 
where 



2Tia' 



62'' 



r) !4 ' 
62 



l + u)h 0(r) -( hm g{T) - g^vD) - - {h{r) - h{r*)) 



g{r) = tan ^(rb) + - log 
h{r) = tan""'^(r6). 



{rby 



((r6)2 + l){rb+iy 



{5.U 



(5.19) 



We can explicitly see that the value of 11^ (r,r) at each point on the string is finite. In 
order to obtain the rate of momentum loss of the quark in this medium, we rewrite (|3.8| ) 
but at this time up to first order: 



1 



-9 



1 



{0)+{l) 



(0) 



(5.20) 



where the subscripts remark the order of perturbation. Some tedious calculation shows 
that the second term in the brackets vanishes on the boundary. So up to viscous order, 
the 1-component of drag force in the RF turns out to be as 



1 



2'Ka' 



62(r) 



1 



1 



1 + b{T) u) 0(r) -h{rl) - -g{rl) 



(5.21) 



One can simply see from F^j^p^ = 7(/^, 0, 0) that the leading contribution in the above 
expression is in complete agreement with ( ^.4|) . 

In analogy with ( p.l4| ), to have a fully covariant expression, one can appropriately 
combine three vectorial objects m^, Mj, 9^0 and construct a general vector which in addition 
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to satisfying (p^), is equal to ( |5.21| ) in the RF. Doing these stages, the result is given by 



27ra' 6^(si 



+ Uf.d(j) I ^h{s2) - ^gis2) - ^ j {{ug.Uf)u'^ + u'}) 

(5.22) 

where Si = x.Uq and S2 = ^—Ug.Uf are lorentz scalars. Scalar nature of the dilaton says 
that demanding to have a plasma with slowly varying coupling in the RF, induces a space- 
time varying dilaton in the LF. The argument of dilaton field in the latter case can be 
determined by the lorentz transformation of the argument in the former case: 



RF ■ 



— > 



LF: (P{-f{t-Ugx)). 



(5.23) 



Such a situation seems somewhat unreal. As a well defined problem, one may expect to be 
encountered with a plasma medium which its coupling varies in time only. We have also 
answered this problem! One can simply observe that in the non-relativistic limit, </> in the 
LF behaves just like in the RF as a function of time. Therefore we can say that drag force 
exerted on a non-relativistic heavy quark moving in a rest viscous fluid, whose coupling is 
varying in time, is given by: 



F\t) 



2Tia' 



62 



km 



Un 



l + ^(|-31og2) 0(t) b{t) 



(5.24) 



In summary, invoking the perturbative procedure explained in § ^ we computed the 
viscous correction to drag force exerted on a non-relativistic quark in a plasma medium 
with a slowly time varying coupling. As it is expectable, by an increase in coupling of the 
plasma, the drag force will increase. 



6 Quark in Bjorken flow 



Bjorken flow [^| has been known as a fine phenomenological description of hydrodynamic 



regime in heavy-ion collisions in last three decades. According to Bjorken's proposal, the 
post-thermalization fluid dynamics in the central rapidity region of heavy-ion collisions is 
longitudinal boost invariant. This symmetry fixes the velocity of the fluid automatically. 
Solving the relativistic hydrodynamic equations of motion determines the fully profile. In 
addition, it is usually assumed that there is no dependence on transverse coordinates. This 
assumption is in correspondence with the large nuclei limit. 

Relevant coordinates to study Bjorken flow are the proper time — rapidity coordinates 
by which the Minkowski metric is be given by 

ds^ = -dr^ + + dx? + dx^ (6.1) 
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a 



Figure 3: A transverse quark in longitudinal expanding plasma 



where r and rj are proper time and space-time rapidity respectively 

1 , t + X3 



The ideal Bjorken flow in this coordinates can be described by 

«'^ = (1,0,0,0), T(r)~-ij. (6.3) 
6.1 Perturbative computations 

At first time, the gravity dual of Bjorken fiow was constructed by JP in the context of 
AdS/CFT. In this derivation, the bulk metric has been computed for asymptotic states, 
namely large times, of the fiow by demanding the Einstein equations to have nonsingular 
solution [jl4| . Studying the energy and momentum loss of a transverse quark in this fiow is 



accounted as a phenomenological problem. This issue has been addressed in |T6|, ^ where 



the authors have similarly applied the JP gravity set-up to compute the drag force at ideal 
order. We choose a different approach to compute the drag force up to viscous order. In 
this way, we exploit the proper time-rapidity coordinates in which, Bjorken flow is at rest. 
As depicted in Fig.|^, we restrict the computation to the case of a quark with zero rapidity. 
For such a moving quark we have 

= 0, T = t. (6.4) 

In order to use our perturbative procedure in § ^, we make a transverse boost to the RF. 
The quark encounters with a global boosted plasma in this frame. Symmetry considerations 
allow us to get the motion in the 1-direction. Consider the quark is moving with velocity 
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Uq in the LF. So in the RF, the Bjorken flow is given by the following profile 



u 



^ = (m°,m\0,0), T{v!'f-u^x^) (6.5) 



where tilded coordinates are related to the RF and = —'jqUq. 

The gravity dual of the above flow can be derived from (|4.1|) as following 



ds^ = —2uq drdr — 2ui da;^dr + (dx^)^ + dr]^ 

(6.6) 



where the derivative corrections are given by 



_ 2 brF{br)u\-ul 



Cii 



3r mqT + Mix^ 
2 uqUi {brF{br) - 1) 
3r uqt + uix^ 
2 brF{br)u1 — u\ 
3r UqT + 



ai = i-^ (6.7) 



In the above expressions, for brevity, we have written b for b{u^f — u^x^). 

In order to construct the string solution in the above gravity background we choose the 
static gauge (cr, r) = (r, f). We assume that 

X^ = e(r,f), X2 = X3 = (6.8) 

where X^(r, f) specifies the embedding of the string in the space-time. The Nambu- 



Goto lagrangian is given by (|3.3| ) and A, B, C are given in [A.4| . Also ^'and ^ are partial 
derivatives of ^ with respect to r and f. 

In here there are some noticeable differences in the drag force computation in compari- 
son within § ^. Firstly, the zero-order embedding of the string is a function of both x^ and 
f through dependency on b. Secondly, since the translational invariance of the action has 
been broken in all boundary directions, the equations of motion in these directions may be 
written as follows 

d£. 

dall- = /„ /, = ^ (6.9) 

where fi are the boundary directions in the RF. As our previous examples, we need to 
evaluate 11^ at the boundary to compute the drag force. So again, we use (|3.4| ) and (|3.9| ) 



by a small difference. We exchange tt^ with 7r(r, r) in this case. We recall from § |^ that 
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at zero order, — AC vanishes somewhere in the bulk, namely the world-sheet horizon. 
Requiring the string to be stretched from the boundary to the horizon implies that the 
denominator should be also vanished at this point. The position of this point will be 
corrected in each order of perturbation. To find the first-order correction, consider 



(6.10) 



where 



(6.11) 



is the zero-order embedding which has been constructed by localizing b in the thermody- 
namic embedding in (|3.12| ). Inserting ( |6.10|) in (|3.4|) , and solving — AC = 0, we obtain 
the first-order corrected position of the world-sheet horizon as r*(f) = TQ(f) +rl{f) where 



r, r 



b 

Uq Ui 



2 b^ r, 



3 U^yf) - j^Uroyf) b'iu'f) 

Tn 



(6.12) 



1 

+ 6 




2 



2uirlF{bTl 





9 1 f Fibr, 
b 



\ ^3 



rlbF{brl) 



UqT 



where the argument of b is u^f in these expressions. As clearly observed in (|6.12|) , rQ(f) is 
the localized version of ( |3.5|) , reconfirmation of our discussion in § ^. It should be noticed 
that the derivative terms have been computed at = 0. Demanding 7r^(r, f) — A to be 
vanished at this corrected r*, the first-order correction to Il\{r,f) at this point can be 
given by 



1 Ml 



27ra' 62 



l + no6(tan-i(r*6)+ ^'^ 



1 + rf 62 



2ni6' tan~^(r* b) 



TT 



71 

2 

2 /'o„,2 



i b (p... _ n2(2n2-l)^ 



(6.13) 



Having n^(r*,f), one can integrate the 1-component of ( |6.9| ) over the r to derive the 
first-order corrected 11^ (r,-?) function. After some calculation we find 



n^(r,f) = nUr*,f)+ / {f,-d,Ul)dT 

tan"^ lb {uib' + upb) + 



^^'^^ '""^ 27ra' 62 



r6 



r262 



{uib' — uob) 



(6.14) 



which is finite at all points on the string. In order to compute the drag force at viscous 
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order, we again utilize ( p.20| ) and obtain 



where 



l+Auob + Buib' + C- 

T 



(6.15) 



^ = 2tan-i(r* b) ■ 
B = tan-\r*o b) + 



2^1 + ro*262' (6.16) 



c = - Kv^+ i^(K)K + Mo)} 

Although the finiteness of corrections is a sign of truth for our perturbative procedure, 
one can separately check it by computing the first-order corrected ^' . After tedious cal- 
culations, it will be turned out that the restriction condition on the bulk embedding of 
the string is satisfied and just, the upper bound of the quark velocity is corrected by a 
first-order contribution. 

So far, we have computed the drag force in the RF. Because of breaking the rotational 
symmetry, we can not covariantize the drag force. So in order to find the drag expression 
in the LF, we first make the inverse boost as following 

four — velocity of the quark : (1,0,0,0) — )■ (7,7^^,0,0) 

1 ^6 1 

quark coordinates : (f , , , 0) — )■ (— t , t , , 0). 

7 

It should be reminded that for a quark with zero rapidity, the proper time and the LF time 
coordinate are coincident. The 1-component of drag force in the LF may be written as 

doc ^ ^ 

For the sake of computing the F^p, we need to find the 11!^. After some calculation similar 
to that for obtaining the 11^, we get that 

F'rf = (6.19) 

up to viscous order. As we have pointed out in § |^, this condition would be maintained 
even when the derivative corrections are taken into account. Physically, this means that 
there should not be any energy loss in the RF at viscous order, just like at ideal order. 
Equation ( |6.19| ) is in complete agreement with this constraint. So substituting (|6.15|) in 



( |6.18| ) gives the rate of momentum loss in the LF. But to proceed further, it is needed to 
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have the relations between the derivatives in these two frames which can be written as 



rf{ 



where u° = 7 and = —'yUq. Finally the rates of energy and momentum loss of the 
in the LF turn out to be as 



dp 
'di 



where 



1 dE 

Uq dt 



A 
B 

V 



-1 Uql 

27ra' &2 



1-6(7^ + 7'^'^) 



(6.20) 
quark 

(6.21) 



2tan-i(y^) -vr, 

vr 77 



tan ^(v/7) 



2 + 7^M 



7 

2„,2 



(6.22) 



+ 7 



1/4 



3 V • ■ 7 
We can write (|6.21| ) in a more clear expression as following 



dp 1 dE Tx 1— 
-dt=-^ = -^^^- 



Uq dt 



Ut 



1 + 



{^'A + i'ulB) 



iiT{t)t 



7D 



(6.23) 



According to our expectation, drag force at ideal order is the same as the drag force in 
thermal plasma medium but with a difference that the temperature would be localized. 
A remarkable point in ( |6.23D is the temperature profile. This is different from the ideal 
Bjorken flow in ( |6.3| ) and should be determined by solving the fluid dynamics equations at 
viscous order. This issue will be discussed in next subsection. 

At the end, we should check the finiteness of the string extension in the 1-direction 
at viscous order. One can do this lengthy and cumbersome calculation and see that ^{r) 
has finite range at viscous order, however, the finiteness of drag corrections is the result of 
allowability of our perturbative method. 



6.2 Phenomenological estimates 

In this subsection we want to apply our perturbative results to the phenomenology of QGP 
in heavy-ion collisions. As it is known, the underlying theory that describes what happens 
through such collisions is QCD. On the other hand, the energy scale in the QGP is related 
to the nonperturbative regime of QCD. Although String/CFT duality could have obtained 
some successes in describing of strongly coupled gauge theories, we need an appropriate 
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dictionary to apply the string/CFT results to QCD. In this way we use the alternative 
scheme introduced by Gubser in |[19|. Using the static force between quarks to normalize 
the 't Hooft coupling, he has compared the strength of the drag force on a heavy quark 
moving through a thermal medium between A/" = 4 SYM gauge theory and QCD. He has 
chosen the energy density (e) as a quantity for comparison. Having esYAi = ^qcd leads to 

TsYM^ ^-^'^Tqcd- (6.24) 



In heavy-ion collisions, the hydrodynamic regime dominates after tth ~ 1/m and the 
expanding hot matter cools down with Bjorken profile. This goes on until the phase transi- 
tion starts at ~ 4/m with ~ 170Mev. During the phase transition, the temperature 
is approximately constant and equal to T^. After the phase transition is over at th ~ 24/m, 
the hadron fluid undergoes a hydrodynamic expansion. Finally at Tf ~ 42 /m the hadron 
breakup, occurs |]26|. But in TV = 4 SYM gauge theory, there is not any phase transi- 



tion. So to get estimate for QCD from this theory, we need to make some assumptions 
about SYM in the range of the temperatures determined by ( |6.24| ). We consider the hy- 
drodynamic regime to be extended from tth to th in SYM. In addition we assume that 
temperature of SYM at th corresponds to Tqcd ~ 170Mev. Applying this assumption to 
( p.'d\) , one can fix the the temperature profile at ideal order 

24 fm 

TsYM{r) = 129Mev (— ^)V3 ij^ < t < 24/m. (6.25) 

T 

The temperature profile at viscous order can be derived from the fluid dynamics equa- 
tions of motion by considering the first-derivative corrections. The width scale of hydro- 
dynamic patches restricts the perturbative computations to the space-time intervals whose 
widths are smaller than 1/T ~ Ifm. Therefore, we use the numerical method to find 
the temperature profile at viscous order. In order to solve the second-order equations nu- 
merically, we fix TsYM = nOMev at r = 24/m. The numerical results has been shown in 
Fig.|. 

Now we are ready to use our perturbative results to obtain estimates on the contribution 
of the drag force on a transverse quark in expanding QGP. To proceed, we consider the 
drag coefficient rjE) which is defined as follows: 

Vd = — -77 (6-26) 
p at 



and to find it in our case, we rewrite the drag force expression in (|6.2ij|) as a function of 
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Figure 4: The numerical profile of the temperature at viscous order. This profile is so 
similar to the analytical profile at ideal order that the deviation could not be simply 
shown. 



the quark momentum: 



dt 2 M ^ ^ 



7rT2(t)'M2 ' Af2^ ' nT{t)V K'P' 



(6.27) 



A, B and V have been introduced in (|6.22|) which should be expressed in terms of p/M 



in the above equation. In contrast to the constant drag coefficient of a quark moving in 
thermal plasma, the drag coefficient is a function of the quark momentum here. Although 
this feature appears only in viscous order, there is no a constant rjo even at ideal order, 
the time-dependency of T leads to variation of the drag coefficient. 

Another important quantity is the diffusion time tD, defined as the time it takes the 
momentum of the quark falls by 1/e. For a charm quark (M=1.5 Gev) moving through 
a thermal plasma with Tsym = l90Mev (and equivalently Tqcd = 250Mev), alternative 
scheme gives tz) ~ 2.1/m [|19]. In our case where the plasma is not in thermal equilibrium, 
one should solve the equation ( |6.27|) to find the ti). To do that, we use the numerical 



profile of the temperature in Fig.|^ and solve the equation ( |6.27|) numerically. In this way 
we need to know the initial time of creating the transverse quark (tq) and whose momentum 
(po) at that time. Respecting our upper-band of velocity in § ^, in Fig.|] we have plotted 
the dependence of the diffusion time on the initial momentum for tq = 1/m case. Just 
like the thermal plasma medium, diffusion time of the quark is independent of the initial 
momentum at ideal order. The viscous curve is below the ideal line, as a direct consequence 
of dissipative effects in a fluid dynamical flow. Another observation in Fig.|^ is that in a 
viscous Bjorken fluid, diffusion time decreases with increasing the initial momentum. This 
is a common behavior in viscous systems. 

In order to complete our discussion, we now study the effect of tq on the diffusion time. 
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Figure 5: Comparison of the diffusion time for a transverse quark created at r = 1/m in 
Bjorken flow between ideal and viscous orders. 
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Figure 6: Diffusion time as a function of the quark creation time, po = 6Gev 



Consider a transverse quark with a high pq near to the upper-band of the momentum. 
Depending on its creation time, it will experience the cooling of the system and the change 
of the temperature. The later creation occurs, the smaller temperatures would be observed. 
The resultant variation of t£, has been shown in Fig.^. Our estimates show that based on 
fluid/gravity duality, the diffusion time of a high energy transverse charm quark could 
be about 4/m, provided that it is created sufficiently late in Bjorken flow. As it can be 
observed, the diffusion time for a charm quark created at tq = 1/m, is about 0.6/m which 
is similar to the result of from the obvious scheme. The origin of dependency on the 
creation time of the quark is that the drag coefficient is time dependent in our case. In 
contrast to previous researches where a thermal plasma was studied, r/z? is a function of 
time in a fluid dynamical flow. 
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7 Discussion 



Applying AdS/CFT correspondence, many problems related to a moving quark in thermal 
plasma have been studied in resent years. The rate of energy and momentum loss 0-0, 0- 
27, jet quenching parameter [^|, Brownian motion [30, 31|, motion in a magnetic 



plasma ||3^ and etc. are all examples of a quark probe which in AdS/CFT duality is 
pictured to a string probe in an asymptotically AdS space-time (for brief review, see |P, |33|). 
In these problems, the solutions of classical string in gravity side, have been strong enough 
to address and clarify some questions in the strongly coupled CFT side. 

It has been accepted that the expanding plasma medium generated after scattering of 
two nuclei evolves from a far-from-equilibrium state to the hydrodynamic regime during 
about Ifm. So one may be interested in studying the above problems but at this time, in 
a fluid dynamical background. To proceed, a string probe should be studied in a gravity 
space-time dual to the boundary fluid dynamical flow. 

In this paper, we have shown that based on fluid/gravity correspondence, one can con- 
struct a perturbative procedure to find the string solution on a perturbed AdS background 
dual to a fluid dynamical flow on the boundary of AdS. According to our discussion in 
§ ^, as far as the string probe lies in one tube in the bulk, the position of the world-sheet 
horizon can be determined perturbatively and so the drag force exerted on the quark will 
take correction, order by order in a boundary derivative expansion. Obviously, by going to 
higher orders in perturbation, the world-sheet horizon position and thereby the drag force 
get higher-derivative correction terms. Requesting the string to be always in a unique tube 
leads to appearance of an upper bound for the relative velocity of the quark and fluid. 
This upper bound is not much lower than unity and is about 0.97. As a remarkable point, 
the necessary condition of being string in its original tube should be checked at each order 
of perturbation, however the finiteness of the computed corrections is a sign of correctness 
in the given order of perturbation. 

By this construction one can study a quark in the fluid medium for a wide range of 
velocities. As an important example we have computed first-derivative correction, namely 
viscous correction, to the drag force exerted on a transverse quark in a QGP experiment. 
An interesting aspect of this procedure is that the drag force at ideal order is exactly the 
localized version of the drag force computed in a global thermal plasma. In another words, 
to find drag force at ideal order, we should insert the local velocity and temperature of the 
fluid in the quark position into ( p.l4| ) at each instance of time. 

On the other hand, existence of the velocity upper bound has a physical outcome to 
our problem. Consider a moving quark with a velocity even grater than the upper bound 
velocity given in ( |4.8| ). Destruction of our perturbative procedure is the first consequence 
of such motion and in addition, we conclude that even at ideal order, the drag force exerted 
on such quark is not a local quantity. It is reasonably related to the high velocity of the 
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quark. We leave further investigation on this issue to the future. 

From the phenomenological point of view, our perturbative method makes improvement 
in previous estimates of quark diffusion time. Our estimates show that in contrast to a 
thermal plasma, the drag coefficient changes to be a time dependent function in a dynamical 
flow. Accordingly, the diffusion time will be a function of the quark creation time. 

This method can open a window to some new problems. Study of rotating quark in a 
fluid or investigating Brownian motion of quark in a fluid dynamical flow are now possible 
through this method. One can generalize |^ to the context of megnetohydrodynamic 
perturbatively, using the metric given in |Q and applying our procedure. Motivated by 
|p7|| , another instructive example is the computation of non-relativistic correction to drag 
force in a fluid medium which can be performed by flnding the string solution in the gravity 
setup of [^. As a complete generalization of our work we leave the study of a quark in a 
general fluid dynamical flow to our further work. 
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A Appendices 

A.l Global boosted plasma 

In this appendix we rewrite the induced metric on the world-sheet given in with more 
detail. From (^) we have 

g = -{Ae + '2Bi' + C) (A.l) 
in which, the coefficients constructed from metric ( p.lD are given by 

A= Gl- GuGu = fipr) 

B = GtrGn - GuGri + u, {GrtGn - GriGn) = 7 {u, - (3) (A.2) 

C = {Grt - U, Grlf = 7' (1 - /3 U^f 

where Uq and /3 are the three-velocity of quark and plasma in the LF respectively. 
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A. 2 Fluid dynamical correction 



Our perturbative computations need to have first-order corrected dual gravity. So in 



following, we recall the computed corrections in |1I5|, ^ up to first order. 



A. 2.1 Free fluid 

The first-order correction to (|4.1| ) has been given in equation (4.24) of [|l^ as follows: 



d.L = rHF(ir) + - r u,u^ a,u>- - r u^d, d.^d.^ ( A.3) 



where 



and 



F(r) 



In 



2 arctan(r) + vr 



(A.4) 
(A.5) 



A. 2. 2 forced fluid 



In this case, the correction to the metric has been given in equation (3.4) of which is 
similar to the ( |A.3p with this difference that the derivatives are covariant derivatives. The 
first-order correction to the dilaton field is given by 



°° - 1 

<l>cor = h u.d(f) I dx 

Irb ^ 



(A.6) 



A.3 Varying plasma 

In this appendix we rewrite all coefficient introduced in ( |5.9| ) with their expressions. Note 
that in following, we have applied first-order corrected metric in 



Ai — G^i — GiiGr 



f{br) 



B\ — G^-t-Gij- — GfiG-j-T — — '"^ 



i?2 — GrrGii — GrlGrl — ^/ 



0^2 



(A.7) 
(A.8) 



and 



C^ = Gl = (4)2 

G2 — GrrG^l — — ^/^/ 

Cs = Gl, = {u)f. 



(A.9) 
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A. 4 Bjorken flow 

The coefficient introduced below (|6.8| ) have the same expressions as which given in previous 
appendix by this difference that inhere our r should be substituted by f. 

References 

[1] J. M. Maldacena, "The Large N limit of superconformal field theories and 

supergravity," Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 
(1999)] [hep-th/9711200]. 

[2] E. Witten, 'Anti-de Sitter space and holography," Adv. Theor. Math. Phys. 2, 253 
(1998) [hep-th/9802150]. 

[3] V. E. Hubeny and M. Rangamani, 'A Holographic view on physics out of 
equilibrium," Adv. High Energy Phys. 2010, 297916 (2010) [arXiv: 1006.3675 
[hep-th]]. 

[4] P. Kovtun, D. T. Son and A. O. Starinets, "Viscosity in strongly interacting 

quantum field theories from black hole physics," Phys. Rev. Lett. 94, 111601 (2005) 
[hep-th/0405231]. 

[5] C. P. Herzog, A. Karch, P. Kovtun, C. Kozcaz and L. G. Yaffe, "Energy loss of a 
heavy quark moving through N=4 supersymmetric Yang-Mills plasma," JHEP 0607, 
013 (2006) [hep-th/0605158]. 

[6] S. S. Gubser, "Drag force in AdS/CFT," Phys. Rev. D 74, 126005 (2006) 
[hep-th/0605182]. 

[7] J. Casalderrey-Solana, D. Fernandez and D. Mateos, "A New Mechanism of Quark 
Energy Loss," Phys. Rev. Lett. 104, 172301 (2010) [arXiv:0912.3717 [hep-ph]]. 

[8] K. B. Fadafan, H. Liu, K. Rajagopal and U. A. Wiedemann, "Stirring Strongly 
Coupled Plasma," Eur. Phys. J. C 61, 553 (2009) [arXiv:0809.2869 [hep-ph]]. 

[9] M. Chernicoff and A. Guijosa, "Acceleration, Energy Loss and Screening in 
Strongly-Coupled Gauge Theories," JHEP 0806, 005 (2008) [arXiv:0803.3070 
[hep-th]]. 

[10] P. F. Kolb and U. W. Heinz, "Hydrodynamic description of ultrarelativistic heavy 
ion collisions," In *Hwa, R.C. (ed.) et al.: Quark gluon plasma* 634-714 
[nucl-th/0305084]. 



30 



N. Banerjee and S. Dutta, "Holographic hydrodynamics: models and methods," 
arXiv:1112.5345 [hep-th]. 

G. Policastro, D. T. Son and A. O. Starinets, "From AdS / CFT correspondence to 
hydrodynamics," JHEP 0209, 043 (2002) [hep-th/0205052]. 

G. Policastro, D. T. Son and A. O. Starinets, "From AdS / CFT correspondence to 
hydrodynamics. 2. Sound waves," JHEP 0212, 054 (2002) [hep-th/0210220]. 

R. A. Janik and R. B. Peschanski, "Asymptotic perfect fluid dynamics as a 
consequence of Ads/CFT," Phys. Rev. D 73, 045013 (2006) [hep-th/0512162]. 

S. Bhattacharyya, V. E. Hubeny, S. Minwalla and M. Rangamani, "Nonlinear Fluid 
Dynamics from Gravity," JHEP 0802, 045 (2008) [arXiv:0712.2456 [hep-th]]. 

G. C. Giecold, "Heavy quark in an expanding plasma in AdS/CFT," JHEP 0906, 
002 (2009) [arXiv:0904.1874 [hep-th]]. 

A. Stoffers and I. Zahed, "Holographic Jets in an Expanding Plasma," 
arXiv:1110.2943 [hep-th]. 

M. P. Heller, R. A. Janik and R. Peschanski, "Hydrodynamic Flow of the 
Quark-Gluon Plasma and Gauge/Gravity Correspondence," Acta Phys. Polon. B 39, 
3183 (2008) [arXiv:0811.3113 [hep-th]]. 

S. S. Gubser, "Comparing the drag force on heavy quarks in N=4 super- Yang-Mills 
theory and QCD," Phys. Rev. D 76, 126003 (2007) [hep-th/0611272]. 

H. van Hees and R. Rapp, Phys. Rev. C 71, 034907 (2005) [arXiv:nucl-th/0412015]. 

M. Rangamani, "Gravity and Hydrodynamics: Lectures on the fluid-gravity 
correspondence," Class. Quant. Grav. 26, 224003 (2009) [arXiv:0905.4352 [hep-th]]. 

V. E. Hubeny, S. Minwalla and M. Rangamani, "The fluid/gravity correspondence," 
arXiv: 1107.5780 [hep-th]. 

S. Bhattacharyya, R. Loganayagam, S. Minwalla, S. Nampuri, S. P. Trivedi and 
S. R. Wadia, "Forced Fluid Dynamics from Gravity," JHEP 0902, 018 (2009) 
[arXiv:0806.0006 [hep-th]]. 

F. Dominguez, C. Marquet, A. H. Mueller, B. Wu and B. -W. Xiao, "Comparing 
energy loss and p-perpendicular - broadening in perturbative QCD with strong 
coupling N = 4 SYM theory," Nucl. Phys. A 811, 197 (2008) [arXiv:0803.3234 
[nucl-th]]. 



31 



J. D. Bjorken, "Highly Relativistic Nucleus-Nucleus Collisions: The Central 
Rapidity Region," Phys. Rev. D 27, 140 (1983). 

K. Yagi, T. Hatsuda and Y. Miake, Camb. Monogr. Part. Phys. Nucl. Phys. Cosmol. 
23, 1 (2005). 

A. Akhavan, M. Alishahiha, A. Davody and A. Vahedi, "Non-relativistic CFT and 
Semi-classical Strings," JHEP 0903, 053 (2009) [arXiv:0811.3067 [hep-th]]. 

M. Ali-Akbari and U. Gursoy, "Rotating strings and energy loss in non-conformal 
holography," arXiv:1110.5881 [hep-th]. 

H. Liu, K. Rajagopal and U. A. Wiedemann, "Calculating the jet quenching 
parameter from AdS/CFT," Phys. Rev. Lett. 97, 182301 (2006) [hep-ph/0605178]. 

J. de Boer, V. E. Hubeny, M. Rangamani and M. Shigemori, "Brownian motion in 
AdS/CFT," JHEP 0907, 094 (2009) [arXiv:0812.5112 [hep-th]]. 

H. Ebrahim and M. Headrick, "Instantaneous Thermalization in Holographic 
Plasmas," arXiv: 101 0.5443 [hep-th]. 

E. Kiritsis and G. Pavlopoulos, "Heavy quarks in a magnetic field," arXiv:1111.0314 
[hep-th] . 

Y. V. Kovchegov, "AdS/CFT applications to relativistic heavy ion collisions: a brief 
review," arXiv:1112.5403 [hep-ph]. 

N. Banerjee, J. Bhattacharya, S. Bhattacharyya, S. Dutta, R. Loganayagam and 
P. Surowka, "Hydrodynamics from charged black branes," JHEP 1101, 094 (2011) 
[arXiv:0809.2596 [hep-th]]. 

S. Bhattacharyya, S. Minwalla and S. R. Wadia, "The Incompressible 
Non-Relativistic Navier-Stokes Equation from Gravity," JHEP 0908, 059 (2009) 
[arXiv:0810.1545 [hep-th]]. 

S. Bhattacharyya, V. E. Hubeny, R. Loganayagam, G. Mandal, S. Minwalla, 

T. Morita, M. Rangamani and H. S. Reall, "Local Fluid Dynamical Entropy from 

Gravity," JHEP 0806, 055 (2008) [arXiv:0803.2526 [hep-th]]. 



32 



